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Abstract. We have investigated the role of transition metal impurities and oxygen 
vacancies in the formation of ferromagnetism in Co-doped Ti02 using LSDA+U 
approach which takes into account strong on-cite Coulomb correlations for electronic 
structure calculations. Several model systems such as supercells of rutile Ti02 
with Co+ 2 ion in high-spin state substituted into titanium cite and with oxygen 
vacancies were considered. We found that exchange interaction of Co magnetic ions 
is ferromagnetic, but very weak due to the large average impurity-impurity distance. 
However, its strength becomes three times stronger when there is a magnetic vacancy 
nearby. The magnetic moment values are 3 \ib for Co 2+ and 1 hb for vacancy with 
the opposite directions of them. Our investigation showed that exchange interaction 
energy of Co and vacancy moments varies from 330 meV to 40 meV depending on the 
distance between them in the supercell. We suppose that strong interaction between 
Co and vacancy moments should be taken into account for the explanation of high 
Curie temperature value in Co-doped Ti02- 
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1. Introduction 

Diluted magnetic semiconductors (DMS's) have been studied extensively because of the 
potential possibility to use both charge and spin degrees of freedom of carriers in the 
electronic devices, namely in the spintronics PHE]. However, the most investigated II- VI 
and III-V compounds doped with magnetic transition metal ions have Curie temperature 
(Tc) about 110 K or less. 

Recent research effort has been focused on developing new ferromagnetic 
semiconductors operating at room temperature. In this connection oxide semiconductors 
became a subject of intent interest El Ej- In 2001 Matsumoto et al. [5, reported 
that Co-doped anatase Ti02 can keep ferromagnetic order up to 400 K with magnetic 
moment about 0.32 /z#/Co for Co concentration up to 8 %. More recently, Park et 
al. [B| have successfully grown the ferromagnetic Co-doped rutile Ti02 films. The 
Curie temperature was estimated to be above 400 K for 12 % of Co content. By now 
a lot of other semiconductors showing room-temperature ferromagnetism have been 
fabricated 0E1E|. 

Despite recent experimental success, there is no consistent theoretical description 
of DMS's. At first, ferromagnetism in Co-doped Ti02 was explained in terms of carrier 
induced mechanism like in III-V based DMS's jH]. Samples fabricated under oxygen- 
rich conditions show negligible magnetization and are insulators [HI E3] so that one can 
suppose that oxygen vacancies, which are easily formed in thin films, are essential to 
ferromagnetic order and conductivity. The magnetic moment value also depends on 
the sample growth conditions. For example, it is reported to be 1.26 /ie/Co in JT] 
to compare with 0.32 /i^/Co, obtained by Matsumoto et al. 5J. The most recent 
measurements ^2] yield a spontaneous magnetization value of 1.1 hb/Co atom for the 
films which were highly insulating with Co in 2+ state. 

The problem of the role of magnetic ions doping and oxygen deficiency in the origin 
of room-temperature ferromagnetism was considered in many theoretical papers. For 
example, in first-principles calculations were performed to consider the influence 

of vacancy position on magnetic properties of Co-doped Ti0 2 , in particular, on Co spin 
state. Park et al. [T3| treated an oxygen vacancy in the supercell (Ti 15 Coi0 3 i) by the 
LSDA+U method with spin-orbit coupling taken in spin-diagonal form. They showed 
that for the anatase structure of Ti 1 _ a .Co x 02 (x = 0.0625 and 0.125) the presence of 
oxygen vacancy near Co results in the state with the spin magnetic moment of 2.53 
whereas the presence of vacancy near Ti does not affect the magnetic moment value 
distinctly. However, according to their results, oxygen vacancy near Ti site is more 
stable. This conflicts with ^1] where the opposite is claimed: the oxygen vacancy 
prefers to stay near Co, and it causes the magnetic moment on Co equal to 0.90 
In jl] overlap polaron model caused by shallow donor electrons was proposed to explain 
high Curie temperatures in Co-doped Ti0 2 . Unlike super-exchange or double-exchange 
interaction, ferromagnetic exchange is suggested to produce long-range order. The other 
recent theoretical works are devoted to the role of interstitial Co that might appear in the 
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material [To"! ITo] . From one point of view, interstitial Co destroys the spin polarization 
of the substitutional Co nearby and, hence, reduces the average impurity magnetic 
moment ^H]- From the other results, there is an enhancement of the average value of 
the local magnetic moment relative to the Co 2+ low spin state with the presence of 
substitutional Co [TH]. Some theoretical investigations claimed that Co impurities have 
strong tendency for clusterization E] • Process of preparation and concentrations of 
Co impurity can lead to Co metal clusters and/or dispersed matrix-incorporated Co [18J. 

At the same time, the authors of El E] confirm the absence of Co metal clusters. 
Thus, there are many aspects of physical properties of Co-doped Ti02 which are under 
debate up to date. 

The main purpose of this work was to clarify the role of Co impurities and oxygen 
deficiency in the origin of ferromagnetism with high Tc in Co-doped Ti02- Since 
tetravalent Ti 4+ ion is substituted by divalent Co 2+ , there are intrinsic oxygen vacancies 
in doped system to compensate this charge difference. The actual value of oxygen 
deficiency is difficult to estimate because defect concentration strongly depends on the 
growth conditions. We have considered several model systems with various content of 
Co impurity ions and vacancies. In our study we assumed that there is no Co clustering, 
Co ions are substitutional and their valency state is +2. Our investigations were 
performed using the linearized muffin-tin orbital (LMTO) band structure calculation 
method within the local spin density approximation with the on-site Coulomb repulsion 
U (LSDA+U) [HHm. The lattice relaxation was neglected. 

The analysis was started with calculation of the electronic structure of Ti02 
in stoichiometric case without doping and defects. Ti02 has three kinds of crystal 
structures: rutile, anatase and brookite. For our calculations we choose the first one 
with crystallographic parameters, taken from [21] . The radii of MT spheres for LMTO 
calculations were Rxi=2.42 a.u. and Ro=1.85 a.u.; two types of empty spheres (atomic 
spheres with zero nuclear charge) were added. The valence band is formed by mainly 
oxygen 2p orbitals contribution, and the conduction band by Ti 3d orbitals. The LDA 
calculation gave the value of band gap equal to 1.7 eV which is significantly less than the 
experimental value of 3.1 eV [3]. Mo et al. [22] shown that while LDA underestimates 
the energy gap, other ground state properties such as equilibrium lattice constants and 
bulk moduli are in good agreement with experimental data. It is well-known problem 
of band gap underestimation in LDA approach, where potential is orbital-independent. 
It can be improved by means of LSDA+U method, where orbital-dependent potential 
acts in different way on the occupied and unoccupied d-orbitals [211] ■ in the top panel 
of figure H the electronic structure obtained from the LSDA+U calculation is shown. 
The values of on-site Coulomb interaction parameters U and Jh for Ti 3c? orbitals were 
chosen to be 6.0 eV and 0.7 eV, respectively [20]. The calculated band gap value is 
3.3 eV in good agreement with the experimental value. The LMTO radii and U and Jh 
values were taken to provide experimental gap value, they weakly influence the magnetic 
properties. 

Below we present results of our LSDA+U calculations for model systems 
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Figure 1. Total densities of states obtained from the LSDA+U calculation for pure 
TiC>2 (top panel) and model oxygen deficient Ti0 2 _i/4 (bottom panel) compounds. 
The figures were aligned by the edge of oxygen bands. Zero of energy is at the top of 
the vacancy band for Ti0 2 -i/4- 



corresponding to 1) Co impurities in Ti02 without oxygen deficiency (section EJ), 2) 
oxygen deficient Ti0 2 but without Co impurities (section EJ) and 3) the simultaneous 
presence of Co impurities and oxygen vacancies in Ti0 2 (section EJ). 

The oxygen vacancies play dual role in this material. In the first, every vacancy 
adds two electrons to the system and in the second, it produces localized states in the 
energy gap. In order to separate the influence of those factors on the electronic structure 
and magnetic state of the system we imitated the first one in the calculation by adding 
electrons to the supercell with compensating their negative charge by the uniform charge 
of the opposite sign distributed over the supercell. This has a meaning of taking into 
account the presence of the vacancy in the crystal but somewhere far away from the 
atoms considered in the model calculations so that one can neglect the influence of the 
localized states produced by this vacancy. 

The substitution of tetravalent Ti 4+ ion by divalent Co 2+ has an effect of adding 
two holes to the system. This could be also imitated in our calculations by removing two 
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electrons from the supercell with the adding the corresponding uniform charge of the 
opposite sign. This adding (removing) of the electrons is not simply shift of Fermi energy 
level as in the "rigid band" approximation because we performed fully self-consistent 
calculations with the changed number of electrons and took into account the change of 
the charge density distribution and hence electron potential. 

In the present study we adjusted total number of electrons in the system so that 
Co ions were always in divalent state (experimental works report El EZH] that all Co 
ions are in +2 formal oxidation state). In each model calculation with Co impurity 
the supercell of Ti02 was considered to make impurity concentration close to the 
experimental values of 8-12 %. For vacancy states we have considered all three possible 
configurations: unoccupied, occupied with one electron and fully occupied with two 
electrons. 

2. Cobalt-doped titanium dioxide 

To investigate the electronic structure of Co-doped Ti02, we used a supercell containing 
sixteen Ti02 formula units with two Ti replaced by Co where impurities are separated 
by 6.59 A. Such supercell corresponds to 12.5 % of Co content. 

In order to reproduce Co ions in +2 formal oxidation state Co 2+ (according to 
experimental evidence 01111123) we have added to the system two electrons per Co ion 
(four electrons per supercell). Physically that is equivalent to considering the system 
with one oxygen vacancy per Co ion (Coi/ 8 Ti 1 „ 1 / 8 2 -i/ 8 instead of Coi/ 8 Ti 1 _ 1 / 8 2 
model compound) but without explicitly taking into account the localized states 
produced by the vacancy in the energy gap. 

The resulting density of states has an insulating character. In figureEJwe present the 
electronic structure for Coi/ 8 Ti 1 _ 1 / 8 2 model compound obtained from the LSDA+U 
band calculations. The values of on-site Coulomb interaction parameters U and J# for 
Co 3d orbitals were chosen as 8 eV and 1 eV, respectively [20] • 

Because of large distances between impurities due to small Co concentration, one 
can expect that exchange interaction between Co ions is weak. To estimate the exchange 
interaction parameter, we examined the total energy difference between ferromagnetic 
(FM) and antiferromagnetic (AFM) orientations of Co ions magnetic moments. Our 
results showed that FM phase is lower in energy than AFM by 3 meV (or ~30 K; that is 
an order of magnitude smaller then experimental value of Tc)- The obtained magnetic 
moments values for Co ions in the ground state are 3 \Lb- Thus, exchange interaction 
of magnetic ions is FM, but weak even for the impurity concentration larger than ones 
in fabricated samples of Co-doped Ti0 2 jEj- Therefore exchange interaction between 
impurities magnetic moments alone without taking into account other states cannot 
yield as high T c as observed. 
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Figure 2. Total and partial densities of states obtained from the LSDA+U calculations 
for antiferromagnetic (AFM; top panel) and ferromagnetic (FM; bottom panel) 
configurations of Co ions in Coi/gTii-i/gC^ model compound. 



3. Titanium dioxide with oxygen vacancies 

In order to explore the electronic structure and magnetic properties changes caused by 
the localized states produced by oxygen vacancies, we at first considered a supercell 
containing four formula units of Ti02 with one oxygen atom replaced by an empty 
sphere; it corresponds to Ti0 2 -i/s composition. Result of the LSDA+U calculation 
is shown in figure fjj (bottom panel). Electronic structure for oxygen vacancy exhibits 
narrow defect band, which is fully occupied. It is formed of mainly Ti rf-states of the 
Ti ions nearby the vacancy split off from the bottom of conduction band and has a 
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Figure 3. Total density of states for oxygen deficient Ti0 2 _i/g model compound in 
magnetic state. Fermi energy is zero. 

symmetry of s orbital centered on the vacancy site. Two electrons added to the system 
as the result of removing O -2 ion occupy this band of s symmetry and result in a 
nonmagnetic ground state. To obtain a magnetic solution and, consequently, localized 
moment, this band should be partially filled. That corresponds to a so-called F + -center 
with one electron localized on an oxygen vacancy state. 

To study half-filled vacancy band (still without explicitly considering cobalt 
impurity), we added one hole per supercell. Physically such occupancy of the defect 
states could be realized if the system contains one Co 2+ impurity per two oxygen 
vacancies or one Ti vacancy per four oxygen vacancies. In all our further calculations 
we added the number of holes which would provide half-filled band formed of oxygen 
vacancies. 

The partial occupancy of the defect band results in a spin-polarized magnetic 
solution. The magnetic state has lower total energy in comparison with nonmagnetic 
one; thus, it is the ground state of the system with the magnetic moment value per 
supercell of 1 \ib- In figure 01 the density of states for this magnetic solution is presented. 

Considering our result one can suppose that not only Co ions, but also vacancies 
yield sizeable magnetic signal in Co-doped Ti02- The similar idea was proposed in recent 
paper [21] , where observed room-temperature ferromagnetism in ZnO was interpreted in 
terms of a spin-split donor impurity-band model, and F + centers, associated with oxygen 
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vacancies, were considered as a possible source of magnetism. In principle, vacancies 
themselves without magnetic impurities could induce magnetization as was shown in 
the work of Monnier et al. |25] . They observed high temperature ferromagnetism in 
normally pure CaB 6 , fabricated under Ca-poor growth conditions. In these crystals 
strong covalency of the B-B bonds makes B-vacancies more energetically preferable to 
form. Large values of magnetic moments could be obtained in the distorted lattices 
where vacancies are located nearby. 

We have also investigated exchange interaction between vacancies, considered FM 
and AFM configurations of vacancies moments. We dealt with a supercell containing 
eight formula units with two oxygen vacancies in a supercell; it corresponds again to 
Ti02-i/s composition. The results of the LSDA+U calculations showed that AFM 
solution is lower in energy by 50 meV. That can be explained by the Hubbard model 
theory. It predicts that AFM solution is the ground state in the case of half-filled band. 

The next step is to study how vacancy moment interacts with Co one and how the 
presence of vacancy influences on Co-Co interaction. 

4. Cobalt-doped titanium dioxide with oxygen vacancies 

In order to explore the influence of vacancy states on the magnetism of Co-doped Ti02, 
we considered a supercell containing eight formula units with one Co impurity and 
one oxygen vacancy at three different possible distances from each other. The energy 
difference between FM and AFM configurations of Co and vacancy was calculated for 
each case within LSDA+U approach to estimate the strength of exchange interaction 
between Co and vacancy magnetic moments. As in the previous section, total number 
of electrons in the supercell was adjusted to produce divalent Co 2+ ions and half-filled 
oxygen vacancy states (F + -centers). That required adding one electron to the supercell 
because every vacancy adds two electrons to the system and two electrons are removed 
when tetravalent Ti 4+ ion is substituted by divalent Co 2+ . That could be realized by 
the presence of additional oxygen vacancy per Co-vacancy pair but without considering 
explicitly the corresponding localized states. 

The obtained results are presented in figure EJ In all three cases, AFM arrangement 
of Co and vacancy spins is lower in energy. The calculated moments of Co and vacancy 
are 3 \ib and -1 /is, correspondingly. Thus, oxygen vacancy prefers to have the magnetic 
moment that is AFM ordered to Co moment. One can see (figure HJ) that strength of 
interaction between Co and vacancy drops rapidly with the increase of the distance 
between them. For the nearest and the longest distances of 1.9 A and 4.1 A, the 
energy difference values are 330 meV and 40 meV, correspondingly. This energy interval 
corresponds to the temperature range from ~3800 K to 500 K. We suppose that the 
presence of strong exchange interaction between Co and vacancy moments might results 
in a strongly coupled system of Co and vacancies magnetic moments with a high ordering 
temperature. 

To investigate vacancy influence on the strength of exchange interaction between 
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Figure 4. The energy difference between FM and AFM configurations of Co and 
vacancy moments as a function of distance between them for Co 1 / 8 Ti 1 _ 1 / s 2 _i/8 
model compound. Distance is given in the units of lattice constant a— 5.92 A. 



Co spins, we used the same supercell as in section |2J but with one oxygen ion removed 
from it. We chose the position of the oxygen vacancy which is at the equal distance of 
3.58 A from both Co ions. Again total number of electrons in the supercell was adjusted 
to produce divalent Co 2+ ions and half- filled oxygen vacancy states (F + -centers). In 
this case that required adding three electron to the supercell. (Similarly to the previous 
case that could be realized by the presence of three additional oxygen vacancies per 
supercell but without considering explicitly the corresponding localized states.) 

Total and partial densities of states for AFM and FM configurations of Co ions 
obtained from the LSDA+U calculations are shown in figure^ One can see that vacancy 
induced states are situated around the Fermi level between occupied and unoccupied 
Co 3c? states without overlapping them. This electronic structure rather agrees with the 
scheme proposed for Mn impurity in a spin-split donor impurity-band model, |21j where 
there is no overlap of vacancy band with 3d states. But in contrast with that, vacancy 
states are not in the middle between Co majority and minority states, but shifted closer 
to the occupied 3d states, although without overlapping, like in the scheme proposed 
for Ti in [21]. 

From the calculations we found that strength of magnetic ions interaction in 
the presence of vacancy is three times stronger than that without vacancy. The 
corresponding energy difference between AFM and FM configurations of Co ions is 
equal to 10 meV; FM orientation of Co magnetic moments is lower in energy. The 
obtained total cell magnetization is 5 //£. At that each Co accounts for 3 /xg. Vacancy 
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Figure 5. Total and partial densities of states for AFM (top panel) and FM (bottom 
panel) configurations of Co ions in Coi/8Ti 1 _ 1 / 8 02_i/i6 model compound obtained 
from the LSDA+U calculations. Zero of energy is at the Fermi energy. 
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Table 1. The list of supercells used in calculations: model of the calculation, type of 
magnetic ordering between Co ions, oxygen vacancies (Vac) or Co-vacancy, A E is the 
difference between FM and AFM configurations. In the third and fourth models the 
distance (D) between Co and vacancy is given. Last two cases correspond to half-filled 
and empty /full- filled vacancy states. 

Model Type of ordering A E, meV 

Co 1/8 Ti 1 _ 1/8 2 _ 1/8 (16 f.u., 2 Co) Co-Co 3 

Ti0 2 _i/ 8 ( 8 f.u., 2 Vac) Vac- Vac 50 

Coi/gTii^/gOa-i/s ( 8 f.u., 1 Co, 1 Vac) Co- Vac (D = 1.9 A) 330 

Coi/gTii.i/gOa-i/g ( 8 f.u., 1 Co, 1 Vac) Co- Vac (D = 4.1 A) 40 

Coi/gTii.i/gOa-i/g (16 f.u., 2 Co, 1 Vac) Co-Co (half-filled) 10 

Coi/gTii^/gCVi/s (16 f.u., 2 Co, 1 Vac) Co-Co (empty/full-filled) 3 



moment is equal to 1 /ie and is AFM ordered with respect to Co moments. The solution 
with FM order of moments is not stable. Again, oxygen vacancy moment prefers to be 
AFM ordered to Co one like in the previous case of one Co and vacancy in the cell, 
which effectively results in a FM coupling between Co ion moments. 

The next step in the investigation of vacancy influence on Co spins interaction is to 
examine it when defect state band is empty or full, and, hence, vacancy is nonmagnetic. 
As vacancy states are situated between Co majority and minority d-states (figure EJ), 
addition of an electron or a hole makes the band full or empty, correspondingly, but 
does not change Co valence. Energy difference between FM and AFM configurations of 
Co ions appears to be about 3 meV in the case of both empty and full defect state band, 
like in the supercell without vacancy. Thus, the presence of oxygen vacancy enhances 
Co ions interaction only when the vacancy is magnetic. 

5. Conclusion 

In summary, using LSDA+U approach we have investigated the role of magnetic 
impurities and oxygen vacancies in the formation of ferromagnetism with high T,^ in Co- 
doped Ti02- We found that exchange interaction of magnetic Co ions is ferromagnetic, 
but very weak due to the large average impurity-impurity distance. However, its strength 
becomes three times stronger when there is a magnetic vacancy present. The obtained 
magnetic moments are 3 /i^ and 1 hb for Co and vacancy, respectively. We also found 
that interaction between Co and vacancy moments is surprisingly very strong even for 
the longest distance of 4.1 A in the supercell. In this case energy difference between 
FM and AFM configurations of Co and vacancy moments is 40 meV, that corresponds 
to ~500 K. Our results lead us to think that the observed magnetic signal cannot be 
attributed to Co ions only: vacancies moments should be taken into consideration too. 
It seems that strong interaction between Co and vacancy moments is a key moment for 
explanation of high Tc in Co-doped Ti02- 
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